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ABSTRACT 

Two  problems  of  memoryless  quantization  and  data  fusion  for  the  detection  of  a  weak 
signal  in  stationary  dependent  noise  are  addressed:  (i)  fusion  from  sensors  with  mutually 
independent  observations  across  sensors  but  dependent  across  time  and  (ii)  fusion  from  sensors 
with  correlated  observations  across  time  and  sensors.  For  each  problem,  we  consider  four  dis¬ 
tinct  schemes  (a)  fusing  the  test  statistics  formed  by  the  sensors  without  previous  quantization, 
(b)  quantizing  suboptimally  each  observation  and  then  fusing,  (c)  quantizing  optimally  each 
observation  and  then  fusing,  and  (d)  quantizing  optimally  each  test  statistic  of  the  sensors  and 
then  fusing  The  observation  sequence  of  each  sensor  consists  of  a  common  weak  signal  dis¬ 
turbed  by  an  additive  stationary  m  -dependent,  <j>-mixing  or  p-mixing  noise  process.  To  guaran¬ 
tee  high-quality  performance,  a  common  large  sample  size  is  employed  by  each  sensor.  Design 
criteria  are  developed  from  the  Neyman-Pearson  test  in  the  fusion  center  for  the  optimal 
memoryless  sensor  test  statistics  and  the  sensor  quantizer  parameters  (quantization  levels  and 
breakpoints);  these  design  criteria  are  shown  to  involve  an  extension  of  the  asymptotic  relative 
efficiency  used  in  single-sensor  detection  and  quantization.  Numerical  results  in  support  of  the 
analysis  are  given  for  the  case  of  dependent  p-mixing  Cauchy  noise. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  N00014-89-J-1375  and  in  part  by  the 
Systems  Research  Center  at  the  University  of  Maryland.  College  Park,  through  the  National  Science  Foundation's  Engineering 
Research  Centers  Program:  NSF  CDR  8803012. 


I.  INTRODUCTION 


Data  fusion  with  multiple  sensors  has  attracted  considerable  attention  in  recent  years  (see 
[l]-[3]),  for  reason  of  cost,  survivability,  and  communication  bandwidth  among  sensors  and  the 
fusion  center.  Most  of  the  existing  results  for  data  fusion  are  based  on  the  assumption  that  the 
observation  sequences  of  the  sensors  are  independent  across  time  and/or  sensors  and  that  each 
sensor  employs  the  scheme  of  a  hard-limiter  type  (threshold  test)  for  its  one-bit  decision.  In 
[4],  the  fusion  scheme  with  a  two-bit  fuzzy  decision  of  each  sensor  is  addressed  for  the  case  of 
a  single  observation  with  dependence  across  sensors.  As  we  know,  in  practice  the  observations 
tend  to  be  dependent  across  time,  when  the  sampling  rates  increase,  and  across  sensors,  if  the 
sensors  are  closely  distributed  geographically.  As  mentioned  in  [4],  the  fusion  scheme  with  one- 
bit  decisions  of  the  sensors  is  too  conservative  and  entails  considerable  information.  We  consider 
optimal  quantization  fusion  in  order  to  utilize  fully  the  channel  capacities  between  sensors  and 
the  fusion  center. 

In  this  paper  we  contribute  to  the  data  fusion  for  dependent  observations  across  time  and 
sensors.  Specifically,  each  observation  consists  of  the  weak  signal  with  a  common  nonrandom 
value  in  additive  noise  with  a  stationary  univariate  density  and  second-order  joint  densities.  In 
addition,  we  compare  the  unquantized  transmission  and  different  schemes  of  quantization  for 
the  transmission  from  each  sensor  to  the  fusion  center.  Two  problems  are  addressed  here:  in 
the  first  problem,  the  noise  processes  of  the  individual  sensors  are  mutually  independent;  in  the 
second  problem,  the  noise  processes  of  the  different  sensors  are  correlated. 

The  model  of  dependence  used  here  is  one  of  the  mixing  sequences,  the  m-dependent,  the 
</>-mixing,  or  the  p-mixing  sequences.  In  general,  for  dependent  noise  across  time  and/or  sensors, 
the  optimal  sensor-fusion  rule  involves  high-order  (larger  than  two)  densities  of  the  noise  that 
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are  difficult  to  characterize,  which  increases  the  complexity  of  the  detection  process.  In  this 
work,  we  adopt  the  suboptimal  scheme  based  on  the  memoryless  nonlinearities  (or  functions 
of  the  observations).  The  scheme  of  optimal  memoryless  nonlinearities  has  been  useful  in  the 
single-sensor  detection  of  a  weak  signal  in  dependent  noise  (see  [5]-[7])  and  has  led  to  the  optimal 
quantizers  in  single-sensor  detection,  if  we  consider  the  quantizers  as  particular  nonlinearities 
(see  [8]  and  [9]).  The  criterion  for  optimal  nonlinearity  or  quantizer  for  each  sensor  is  derived 
from  the  Neyman  Pearson  test  in  the  fusion  center. 

In  both  problems  described  above,  four  schemes  based  on  nonlinearities  and  quantizers  are 
considered  with  a  common  large  sample  size  for  each  sensor,  in  order  to  form  the  test  statistics 
of  the  fusion  center.  Note  that  using  a  large  sample  size  guarantees  a  high-quality  performance 
in  detection  of  a  weak  signal.  In  the  first  scheme,  the  test  statistic  T„^,  which  consists  of  the 
sum  of  n  observations  Xjk\  ( i  =  1,2, •••,»)  passes  through  an  appropriate  nonlinearity  gk(-), 
i.e.,  Tn,k  =  5Zr=i  9k(Xlk)),  is  formed  by  the  k- th  sensor,  and  transmitted  directly  to  the  fusion 
center,  where  a  likelihood  ratio  test  is  performed.  In  the  second  scheme,  the  observation  Xj^  of 
the  A-th  sensor  at  time  i  is  first  quantized  by  a  quantizer  gk(-)  obtained  from  the  discrete-form 
of  the  optimal  nonlinearity  gk  of  the  first  scheme  and  then  transmitted  to  the  fusion  center, 
where  the  test  statistic  Tn<k  —  H"=i  gk(xj^)  is  formed  by  collecting  the  quantized  data  n  times 
from  each  sensor  and  a  likelihood  ratio  test  is  performed  on  the  basis  of  Tntk,k  =  1,2, 

In  the  third  scheme,  a  scheme  like  the  second  one  is  employed,  except  that  the  quantizers  are 
now  optimized  and  not  related  to  the  nonlinearities  of  the  first  scheme.  In  the  forth  scheme,  the 
test  statistic  of  each  sensor  is  first  formed  as  in  the  first  scheme  by  a  nonlinear  preprocessor  and 
then  quantized  before  transmitting;  in  addition,  we  show  that  the  complex  optimization  for  the 
optimal  nonlinearity  and  quantizer  can  be  separated  in  two  decoupled  optimization  problems, 
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which  considerably  simplifies  the  original  optimization  problem. 

As  mentioned  above,  the  design  criterion  for  the  optimal  nonlinearity  and  quantizer  for  each 
sensor  is  from  the  Neyman  Pearson  test  in  the  fusion  center.  It  turns  out  that  the  final  form  of 
the  design  criterion  is  the  deflection  in  that  the  test  statistics  are  jointly  Gaussian  distributed 
under  a  large  sample  size  (see  [10]). 

The  remainder  of  this  paper  is  organized  as  follows:  in  Section  II,  the  weak  signal  model  is 
described  and  basic  assumptions  are  given;  in  Section  III,  the  case  of  dependence  across  time 
only  is  considered;  in  Section  IV,  the  case  of  dependence  across  time  and  sensors  is  considered; 
in  Section  V,  numerical  results  from  simulations  are  given  to  illustrate  the  analysis;  in  Section 
V,  the  conclusions  are  drawn. 
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II.  PRELIMINARIES 


The  problem  we  address  is  modeled  as  a  distributed  hypothesis  testing  in  a  multi-sensor 
environment  involving  K  sensors 

II  ^ 

Il\k)  :  XiW  =  9  +  N{(k\  i  =  k  =  (2) 

where  9  =  C f  \/R  — *  0,  as  R  —*  oo,  and  C  is  a  constant.  For  each  sensor  k,  {N^k\  i  =  1,  •  •  • ,  n} 
is  a  stationary  dependent  noise  sequence  with  a  univariate  density  fk(-). 

Let  gk  be  a  nonlinear  function  (nonlinearity)  of  X-k^  for  the  sensor  k  and 


fik, o  =  E0[gk(xj%  gk<6  =  Ei[gk{xlk)] 


(3) 


be  the  means  of  the  nonlinearity  for  sensor  k  under  Ho  and  H\ ,  respectively.  Throughout  this 
paper,  we  assume  that  appropriate  conditions  are  satisfied  for  the  mixing  (the  m-dependent,  the 
4>-  or  the  p-mixing)  model  of  dependence,  such  that,  under  Hi  ( i  =  0, 1),  for  sensor  k  and  all 
9  >  0, 

m 

°lfi{9k)  =  t>ar*[0*(-Xifc))]  +  2  £  cove[9k{X[k))gk(X jj\)]  >  0  (4) 

j=i 

converges  absolutely  as  m  —>  oo,  and  under  the  large  sample  size  n, 

A=  YHvk{X(jk))  -  Hk,e{gk) ]  (5) 

V71  j-i 

is  asymptotically  Gaussian  distributed  with  mean  zero  and  variance  a\  6(gk ),  according  to  the 
central  limit  theorem  (see  [11-13]).  We  also  make  the  following  regularity  conditions,  for  9  — *  0 
and  under  a  large  sample  size 

J  9k(x)fk(x  -  6)dx]  =  J  J jj[gk(x)fk(x  -  9)]dx  (6) 
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(?) 


fik,e(9k)  -  Hk,o(9k)  dfie(gk) 

e  oe 

<Tk,e{9k)  -*  °k,o(9k),  (8) 

and  Hk,e  >  (*k, o,  where  the  integration  is  over  the  range  (— oo,  oo);  this  range  is  used  for  integra¬ 
tions  in  this  context,  when  there  are  no  explicit  indications  for  the  ranges  of  integrations. 

Under  the  above  framework,  conditioned  on  if,  (i  =  0,1)  and  for  each  sensor  k,  Tn<k  = 
gk(X-^)  is  asymptotically  Gaussian  distributed  with  mean  ngk,$(9k)  and  variance  na\  0(gk), 
for  all  9  >  0.  Furthermore,  ( TUik ,  k  =  1,2,  •  •  •. K )  are  multivariably  Gaussian  distributed. 

In  the  following  sections,  we  omit  the  argument  g^  of  fik,e{‘)  and  for  all  9  >  0  for 

the  purpose  of  convenience;  the  m  is  finite  for  m-dependent  noise  and  m  — +  oo  for  (^-mixing  or 
/>-mixing  noise. 
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III.  THE  CASE  OF  CORRELATION  ACROSS  TIME 

In  this  section,  we  consider  the  case  of  dependence  across  time  only,  i.e.,  {X,^}”=1  and 
{A^}"=1  are  mutually  independent  for  k  ^  j.  The  four  schemes  mentioned  in  Section  I  are 
addressed  and  optimal  nonlinearities  as  well  as  quantizers  are  formulated. 

III. A.  Unquantized  Transmission 

We  first  consider  the  scheme  described  in  Figure  1  in  which  each  sensor  k  transmits  a  test 
statistic  Tn<k  =  ]C?=i  )  to  the  fusion  center,  where  the  Neyman  Pearson  test  is  used  on 

the  basis  of  the  transmitted  data. 

Let  In  L  be  the  natural  log-likelihood  ratio  function  of  Tntk,  k  =  1,2,**-,  JTsT .  The  test  of  the 
fusion  center  for  the  fixed  level  a  is  described  as  follows 


P0(ln  L  >  r])  =  a 


min  Pi(h\L<ri) 

gk,k=l,2,-,K 


(9) 


where  rj  is  the  threshold  determined  by  a.  Since  Tnik,k  =  1,2, are  independent  for 
different  sensors,  the  likelihood  ratio  function  of  TUtk,  k  =  1, 2,  •  •  •,  K  has  the  following  form 

i  T  _  1  fl<,i(Tn,k ?  k  =  !,•••,#')  _  v-'i  fk,i{Tn,k)  _  ~  Mfc.Orp  ^k,e  ~  Vk, o„i 

III  iy  111  „  rr\  2—J  111  A  (rV  \  Z-A  -2  1Tl1k  0_2 


fK,o(Tn,ky  &  —  1  ?  *  ”  i  -ft  )  fk,o(Tnyk)  j.=1 


whose  expectations  and  variance  under  two  hypotheses  are 


Ho 


and 


£i  [In  L\  =  -£0[ln  L]  =  n  £  ^  ^,o)2 

ti  Ho 


Var[lnL)  =  nJ2—^—2 


(10) 


(11) 


(12) 


fc=l  fc,o 

As  mentioned  in  Section  II,  Tn>fc  for  each  k  is  Gaussian  distributed  under  either  hypothesis 
and  the  large  sample  size  n.  Therefore,  from  the  form  of  (10),  In  L  is  Gaussian  distributed  under 
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either  hypothesis  and  a  large  sample  size,  and 


„  ,,  r  n  ,ln  L  -  .Eh [In  L]  r)  —  £0[ln  L]  N 

p°i]n  L>')= p°(  >  =  rI 


V  -  £0[ln  L] 


sJV  ar  [In  L\ 


=  a 


(13) 


where  erfc(-)  is  the  complement  error  function, 


7]  =  \Jv ar  [In  L\erfc  x(a)  +  Eo[\a  L] 


(14) 


and 


_  „  .  In  L  —  E\  [In  L]  rj  —  E\  [In  L] . 

P,(lnl<  n)  =  fi(  w  <  '/.r-nl-n1) 


x/FarjTnZJ  \/V  ar[ln  L 


=  er/c 
=  er/c 


Ei  [In  L]  -  t) 


=  er/c 


[v/Kar  [In  L\  J 

(/^fc,6>  -  Mm) 


£i[ln  L]  -  £0[ln  L] 


—  erfc  J(a) 


»E 

<  k- 1 


I<  ...  n2\X/2 


^/Varpn  X] 
—  er/c-1 (a) 


'fc,0 


(15) 


With  this  form  of  miss  probability,  the  minimization  problem  (9)  is  characterized  by 

(£i[ln  X]  -  -Epftn  L])2 
max  _  _ 

gk,k=i,2?~,K  Var[  In  L ] 

V-''  (M k,B  ~  fJ-kfi)2 
—  n  ■  max  >  - — 1 — ^ — -L- 

gk,k=l,2,..;K  ^ 


(16) 


k= 1  ''M 

for  fixed  a,  where  (£i[lni]  —  JSofln  L])2/Var[ln  L]  is  the  deflection.  Then  under  the  regularity 
condition  (7), 

(17) 


v~ '  (MM  —  Mfc, o)2  n2  V'  (^k, o)  o 

„2  “*■  ^  2_/  -2  - 


fe=i  aM 


fc=l  fc,0 

Since  0  is  a  constant  and  the  observations  are  independent  across  sensors,  the  optimization 
of  (17)  with  respect  to  g t,  for  =  1,2, •  •  -,K  can  be  conducted  separately,  i.e., 


(Aflf 

alo 


max 


(18) 


for  k  =  1,2, •••,/('.  This  minimization  problem  has  been  solved  for  m-dependent  noise  in  [5] 
and  for  (/-mixing  noise  in  [6]  (see  also  [14]).  Here  we  only  give  the  final  form  of  the  integral 
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equation,  which  determines  the  nonlinearity  gk(-),  for  k  =  1, 2  •  •  • ,  K,  as  follows 

-  f'k(x)/fk(x)  -  J  Kk(x,y)gk{y)dy  =  gk(x),  k  =  1,2, — ,  if  (19) 

where  /£(•)  is  the  derivative  of  /*.(•),  and 

m 

Kk{x ,  y)  =  2  J2  f^ll/Nl(y\x)  -  (2 TO  +  l)/*(y)  (20) 

j=l 

is  the  kernel  of  integration  with  the  conditional  density  /^^^(j/lx). 

III.B.  Directly  Quantized  Transmission 

Under  the  scheme  of  III.A,  each  sensor  has  to  transmit  Tn<k  =  IZiLi  ffk(xj^))  (k  =  1, 2,  •  •  • ,  A'), 
a  real  number,  to  the  fusion  center,  which  is  inpractical  from  a  bandwidth  point  of  view.  The 
easiest  way  to  modify  the  structure  of  III.A  is  to  approximate  the  the  integral  equation  (19), 
which  is  satisfied  by  the  optimal  gk{‘),  with  its  counterpart  of  a  discrete-form  and  let  each  sen¬ 
sor  transmit  the  finite-dimensional  nonlinearity  gk(X^)  (i  =  1, 2,  •  ■  n)  of  the  discrete-form  of 
optimal  gk( •)•  Then  the  test  statistic 

TU  =  (21) 

f=l 

for  each  sensor  is  formed  in  the  fusion  center.  In  fact,  this  modified  structure  turns  the  nonlin¬ 
earity  of  sensor  fc  to  a  quantizer  and  gk  stands  for  the  breakpoints  and  the  quantization  levels 
(discrete- function  values  of  gk  at  breakpoints). 

In  order  to  approximate  the  integral  equation  (19)  by  its  discrete-form,  we  set  the  interval  of 
integration  to  be  (X^]n,  X^lx)  (chosen  according  to  the  main  coverage  of  the  noise  density)  for 
each  sensor  k  and  write  the  discrete-form  of  an  M-level  quantizer  for  each  sensor  k  as  follows 

M 

-  A(4‘,)/A(*i‘))  =  »(4?)  (22) 

i=o 
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where  x£f\i  =  0,1,2 ,  (x£^  =  X^-n,x{k2  =  X$ax)  are  the  breakpoints  of  the  quantizer 

and  Ax£J°  is  the  approximation  of  dxj  for  sensor  k,  (for  example  Axi^  =  (X^L  -  X^}n)[M). 
Define  the  vectors 

st=  (23) 

fl  =  [/«4S>m4fV--,/*(zSi!)]  (24) 

and  the  matrix 

Ok  =  |gS‘>]  (25) 

where 

=  A(zi‘l)(A'i(z‘‘»,x<;»)Aa!l*l  +  S(U  -  (,)].  (26) 

Then  the  algebra  equation  (22)  can  be  written  as 

fiT  =  Gk9kT  (27) 

for  each  sensor  k  (k  -  1,2 ,---,/v).  To  solve  (27)  for  g^,  we  assume  that  for  each  k,  x{k) 

( i  =  0, 1,  •  •  •,  M)  are  chosen  such  that  the  matrix  Gk  is  nonsingular.  Thus  we  have 

gj  =  Gk~'f£  (28) 

for  k  =  1,2 From  the  discrete-form  ^(x^)  (i  =  0, 1, •  •  •, M)  we  can  characterize  the 
quantizer  for  sensor  k  and  any  observation  x  as 

9k{x £?)  if  x  <  xik) 

M*^)  +  gk(x(J:l)}/2  if  x  6  (^,4.]  (29) 

9k{xikM- 1)  if  x  >  x£^]_j 

for  i  =  1,2,  -  2. 
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In  other  words,  each  sensor  k  transmits  the  quantized  value  gk(x\k^)  described  above  n  times 
to  the  fusion  center,  where  the  test  statistic  Tntk  is  formed  and  a  likelihood  ratio  test  (described 
in  Section  III.  A)  is  performed  on  the  basis  of  the  Tn<k ,  for  k  =  1, 2,  •  •  • ,  K . 

Note  that  in  the  above  quantization  scheme,  the  number  of  levels  M  and  the  breakpoints 
a:^  ( i  =  0,1,---,M)  are  not  necessarily  the  same  for  distinct  sensors.  For  the  purpose  of 
simplicity,  we  may  set  them  to  be  the  same  for  all  sensors.  Although  the  resulting  quantization 
scheme  will  not  be  optimal,  its  performance,  as  derived  from  simulations,  is  acceptable  for  a 
reasonable  number  of  quantization  levels.  Moreover,  the  analysis  and  the  design  of  this  scheme 
is  straightforward. 

IIIC.  Optimally  Quantized  Transmission 

In  the  scheme  of  Subsection  III.B,  the  quantizer  ffk  of  sensor  k  is  obtained  directly  from  the 
discrete-form  of  the  optimal  nonlinearity  gk  of  the  first  scheme,  and  thus  is  not  an  optimal  quan¬ 
tizer.  Since  quantizers  also  function  as  nonlinearities  (thus  satisfing  the  central  limit  theorem 
and  regularity  conditions  of  Section  II),  we  can  reformulate  the  deflection  of  each  sensor  for  the 
second  scheme  by  initially  working  on  the  nonlinearity  of  a  quantizer-form  and  by  deriving  the 
associated  deflection  as  the  design  criterion  for  the  optimal  quantizer  (breakpoints  and  levels) 
of  each  sensor.  Actually  this  type  of  scheme  has  been  addressed  in  [8]  and  [9]  for  single-sensor 
detection.  By  using  the  technique  similar  to  the  one  of  [9],  we  can  obtain  the  optimal  quantizer 
for  each  sensor  through  numerical  optimization  methods. 

Let  Qk(X^)  be  the  quantizer  with  M  levels  for  sensor  k  (k  =  1,2,  •  •  •,  K)  in  our  scheme. 
Denote  by  4  =  [4, o  4,1  •  •  •  4,m]  its  breakpoints  and  by  Uk  =  Ufc)2  •  •  • Uk,M }  its  levels. 
The  test  statistic  formed  from  the  transmission  of  sensor  k  now  is  Tn>k  =  HILi  Qk(X-^).  As 
mentioned  in  Section  II,  we  assume  that,  under  the  large  sample  size  n,  each  test  statistic  is 
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Gaussian  distributed  and  the  Neyman-Pearson  criterion  in  the  fusion  center  is  equivalent  to 
maximizing  the  deflection.  In  addition,  since  Qk  of  the  sensor  k  is  one  special  form  of  the 
nonlinearity  described  in  Section  II,  the  assumptions  and  regularity  conditions  are  still  useful 
here. 

0 

For  this  case,  the  deflection  in  the  fusion  center  has  the  form 


K 


K 


Dl<,$  =  ^2  DW  ~  H 

k=l  fc=l 


( njlk,e  -  njlkfi)2 

n*l,o 


(30) 


where  nfik,6  and  na^fi  are  the  mean  and  the  variance  of  Tn^  for  all  Q  >  0,  and  thus  functionals  of 
Qk-  Again  we  omit  the  argument  Qk  in  the  formulations  involving  the  means  and  the  variances. 
From  (4)  we  have 

M  (  m  M  M 

»kfl  =  £(uk,i)2  /  ’  fk(x)dx  +  2j2YUluk,iuk,ipo{Xik)  e  ■X'JJ.I  €  (tk,i-i,tk,i]} 

l=i  j= i  i= i  i=i 

—(2m  +  /  ’  fk(x)dx)2  (31) 

;=i  dtk,  i_i 


The  deflection  criterion  is  the  optimization  problem 


Q>;£%,.,KDKMk)- 


(32) 


Since  the  observations  are  independent  across  sensors,  the  above  optimization  problem  is  decou 
pled  for  different  sensors  and  turns  out  to  maximize 

(Am  “  M2 


°l,o(Qk) 

with  respect  to  Qk ,  as  6  —*  0.  Equivalently,  we  maximize  the  efficacy 

1-  DkAQk)  _  F  (p  x  _  {dfik,e/de)2\g=0 

with  respect  to  Qk  individually,  for  k  =  1,2, ••  -,K-  In  (34), 
dfik,e 


de 


r  M 

■|e=o  =  J  Qk(x)fk(x)dx  =  J2  uk,i[fk{tk,i- i)  -  fk(tk,i)] 


(33) 


(34) 


(35) 
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The  optimization  problem  of  (34)  has  been  addressed  in  [9]  for  the  single-sensor  detection 
with  the  m-dependent  noise  of  symmetric  probability  density;  a  method  similar  to  that  in  [9]  can 
be  used  to  search  for  the  optimal  levels  uk,i  (Z  =  1,2,--  •,  M)  and  breakpoints  tk,i  (l  =  0, 1,  ■  •  •,  M) 
for  noise  with  a  asymmetric  density.  We  do  not  duplicate  the  manipulations  of  [9]  and  give  only 
the  results. 

For  k  =  1, 2,  ■  •  • ,  I{  and  /  =  1, 2,  •  •  • ,  M,  define 

&flk)  =  (36) 


the  vector 


the  matrixes 


with 


and 


with 


AA  =  [A/1WA/f)...A/iJ)] 


Fk  =  diag{  /  fk(x)dx  ■  ■■  fk(x)dx} 
Jt0  dtM-l 

ft  =  Kj1] 


i=i 


Rk  =  [<?] 


1$  =  f  fk(x)dx  •  f  fk{x)da 
dti- 1 

Then  the  efficacy  of  sensor  k  has  the  form 

Ek(Qk )  - - ^  - 

VW  u^(Fk  +  Pk-Rk)uk 


(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 
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In  this  scheme,  we  assume  that  tkri  ( l  =  0, 1,  •  •  •,  M)  of  sensor  k  are  chosen,  such  that  the  matrix 
F  +  P  —  (2m  +  1)-R  is  positive  definite.  Then  the  optimal  quantization  levels  of  sensor  k  are 

uk  =  ~{Fk  +  Pk-  Rky'Afk  (44) 

for  fixed  breakpoints.  Substituting  the  optimal  levels  given  by  (44)  into  the  form  of  Ek,  we  have 
the  efficacy  of  the  sensor  k  as  a  function  of  the  breakpoints  only,  i.e., 

Ek  =  A fkT(Fk  +  Pk-  Rky'Afk  (45) 

where  we  omit  the  arguments  //  (/  =  0,1,--*,M).  Then  numerical  optimization  techniques, 
such  as  the  gradient  method,  can  be  used  to  solve  the  optimization  problem  of  (45)  and  obtain 
first  the  optimal  breakpoints  and  then  the  optimal  levels. 

III.D.  Optimally  Quantized  Transmission  with  Nonlinear  Preprocessors 

The  most  difficult  part  of  the  scheme  in  Subsection  IIIC  is  the  multi-dimensional  optimization 
problem  of  (45),  which  is  time-consuming  in  the  case  of  a  large  m  when  we  are  searching  for 
the  stationary  points,  particularly  when  the  number  of  quantization  levels  M  is  large.  In  this 
section,  we  consider  a  modified  scheme,  in  which  a  nonlinear  preprocessor  is  used  for  each  sensor 
to  form  the  test  statistic  and  handle  the  dependence  before  quantizing.  This  scheme  leads  the 
search  for  optimal  quantizers  to  a  special  case  of  the  optimization  problem  presented  in  the 
previous  section  for  m  =  1,  then  saving  much  time  during  the  searching  procedure.  In  Figure  2, 
we  give  the  configuration  of  this  scheme. 

Each  preprocessor  consists  of  a  continuous  nonlinear  function  fifjt(-)  (k  —  1,2,  •  •  •,  K)  of  its 
observation  X ,  ( i  =  1,2,  •••,«).  The  function  of  the  preprocessor  k  is  the  same  as  in  the  first 
scheme,  namely  to  form  a  test  statistic  {Tn<k)  by  passing  its  observation  through  the  nonlinearity 
gk  and  summing  up  to  the  large  sample  size  n,  as  that  Tn<k  are  Gaussian  distributed  individually. 
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Following  each  preprocessor,  a  quantizer  qk  is  designed  for  the  quantized  transmission  from  the 
sensor  k  to  the  fusion  center.  Then  the  Neyman- Pearson  test  is  used  in  the  fusion  center  based 
on  these  quantized  data.  In  this  scheme,  we  also  assume  a  large  number  of  quantization  levels 
M,  so  that  the  quantized  data  qk(Tnjk )  of  a  Gaussian  random  variable  Tn>*  are  still  Gaussian 
distributed.  The  optimal  nonlinearities  gk  ( k  —  1,2,  •  •  • ,  K)  and  the  companioned  optimal 
quantizers  qk  are  obtained  by  using  the  deflection  criterion  from  the  Neyman-Pearson  test  of 
Gaussian  data. 

Let  Mk,e,  Vk0  be  means  and  variances  of  the  quantizers  qk(Tn,k)i  for  k  =  1,2,  •  •  •,  K,  under 
Hi,  for  9  >  0,  or  H0,  for  9  —  0;  both  of  them  are  functionals  of  ( gk,qk ),  which  are  omitted  in 
the  following  formulations.  Thus, 


Mk,e  =  J  qk(x)fG( )dx  (46) 

and 

Vk,6  =  J  -  ( J  qk(x)fG(X--^-^  )dxf  (47) 

where  /g(-)  is  the  standard  Gaussian  density,  and  nfik,e  and  nak  0  are  the  mean  and  the  variance 
of  THtk,  for  all  9  >  0,  and  thus  functionals  of  gk- 

As  in  the  previous  section,  the  deflection  criterion  is  equivalent  to  maximizing  the  efficacy 
of  each  sensor  k,  as  9  — ►  0,  because  of  independent  observations  across  sensors.  The  efficacy 
of  each  sensor  is  now  a  functional  of  gk  and  qk,  and  the  design  criterion  is  to  maximize  it  with 
respect  to  (gk,qk)>  for  k  =  1,2,  •  •  -  ,I(  and  0  — >  0,  i.e., 


maX9k,qk 


{dMk,e/d9}2\e=o 

v ,2 

vk,  0 


(48) 


In  the  following,  we  show  that  the  optimization  problem  (48)  for  each  sensor  k  can  be  separated 
into  two  decoupled  optimization  problems. 
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Changing  the  variable  and  using  a  condition  similar  to  (6)  we  have 

Mkft  =  J  Qk{x)foC ^0~)dx  ~  J  QkiVnvkfiX  +  nnk,e)fG(x)dx  (49) 

and 

=  J  \Zv.rj t  /yr.  +  nnt,e)la(^)dx 

=  Vn  itM.  f  [  »  +  „,,M)]/0(z)<fo 

<rjt)0  do  J  ox 

=  •  /  qkWwkfi*  +  n»k,e)fG(x)dx  (50) 

^A:,0  •/ 

The  optimization  problem  of  (48)  is  characterized  by 

Mo(gfc)  \2  _ (/ qkiVnOkfiX  +  ni.ikto)f'G(x)dx)2 _ _  (51) 

<Jkfi(§k)  f  qk(\/^<7k,ox  +  n/.ikfi)fG(x)dx  —  (/ qk{'/™ak,ox  +  nfik,o)fG(x)dx)2 

Let  $*(*)  =  qk(okfix  +  W,o),  for  k  =  1,2,  •  •  •,  JRf.  Then  the  optimization  problem  of  (51)  can  be 

conducted  separately  as  follows.  We  first  maximize  the  objective  function 

q _ U  qk(x)fb(x)dxy _  f52) 

qk  I  ql(x)fG{x)dx  -  (f  qk(x)fG(x)dx)2 

with  respect  to  the  quantizer  qk  (i.e.  the  breakpoints  and  the  quantization  levels)  for  fixed  <yk,o 
and  fik,o ,  which  is  a  special  case  of  the  previous  section  for  m  =  1.  Because  of  the  form  of  Sqk, 
the  optimal  breakpoints  and  the  quantization  levels  depend  only  on  /g(-)  and  are  independent 
of  (ok,o,Hk,o),  and  thus  gk.  Then  the  optimal  gk  can  be  obtained  by  maximizing  the  objective 
function  S9k  with  respect  to  gk,  which  has  been  shown  in  Subsection  III. A. 

After  the  optimal  gk  and  qk  are  determined  for  all  sensors  k  =  1,2,  •  •  -,K,  we  can  evaluate 
ak,o(gk),  as  well  as  Hk,o{gk);  the  optimal  quantizers  qk  are  implemented  as  qk(x)  =  qk[(x  - 
l*kfl)l*kfl]*  for  all  x. 
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IV.  THE  CASE  WITH  CORRELATION  ACROSS  TIME  AND  SENSORS 


We  now  consider  the  dependence  across  time  and  sensors  for  the  case  of  two- sensor  detection. 
Although  we  do  not  address  the  general  case,  the  results  of  this  section  can  be  extended  to 
multiple-sensor  detection  with  more  than  two  sensors. 

IV. A.  Unquantized  Transmission 

The  sensor-fusion  scheme  employed  here  is  the  same  as  the  one  in  the  Subsection  III.A,  i.e., 
each  sensor  passes  a  test  static  Tn<k  =  ]C”=i  (k  =  1,2)  to  the  fusion  center,  where  the 

Neyman-Pearson  test  is  performed.  ■< 

Let  po  be  the  correlation  coefficient  of  Tn>\  and  Tn> 2,  under  Hq.  As  discussed  in  Section  II, 
since  the  signal  is  weak,  it  is  reasonable  to  assume  that,  under  either  hypothesis  Hi  ( i  =  0, 1), 
the  correlation  coefficient  of  Tn>  1  and  T„; 2,  for  all  0  >  0  and  0  — >  0,  has  the  same  form 

Eo[(Tnii  —  np,\to)(Tn,2  - 


pe(gi,g  2)  —  poigugi)  = 


(53) 


yJVar(Trh\  )Var(l\tt2) 

which  has  been  studied  in  [15].  In  the  following,  we  omit  the  arguments  of  p  for  greater  conve¬ 
nience.  The  likelihood  ratio  function  of  Tn<k  ( k  —  1,2)  in  the  fusion  center  has  the  form 
fK,l(Tnt\,Tn<2) 


In  Z=ln 
1 


1  -Pi 


/jf,o(rntl,rat3) 

Pi, 8  ~  Pl,0  Po(P2,8  ~  P2,0 ) 


'1,0 


al,  0^2,0 


Tn,\  + 


P2,8  ~  P2,0  Po{Pl,8  ~  P l,o) 


a. 


2,0 


<yXI0°'2,0 


'  n, 2 


A4? ,8  ~  ^1,0  p\,e  ~  P 2,0  ,  po{Pl,9P2,e  -  Pl,oP2,o)-i 

o,2  +  -  T 

z<72,0 


(54) 


2^0  Z<r2,0  <^1, 0^2,0 

It  is  straightforward  to  show  that  the  expectations  and  the  variance  of  In  L  under  two  hypotheses 


are 


E1[\nL\  =  -E0[ln  L] 
1 


2(1  -Pi) 


n{pi,9  -  Pi,o)2  .  n(p2,9  -  P2,o)2 
•>  '  2 


1,0 


'2,0 


2np0(plie  -  pifi)(p2,9  -  P2, 0) 
<71,002,0 


(55) 
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and 


VarM=(WS) 


1  [ n(n a,,?  -  Mi,o)2  n(M2,®  -  M2, o)2  2n/?0(Mi,g  -  Mi,o)(M2 ,®  -  M2,o) 


^1 ,00-2,0 


Following  steps  similar  to  those  in  Subsection  III.A  for  the  Neyman-Person  test  we  obtain 
the  deflection  criterion 


D{L)  = 


(£i[lnX]-£0[lnX])2 
Far  [In  X] 

n  (MM  ~  Mi,o)2  |  (M2 ,®  -  M2,o)2  2ft0(Mi,g  -  Mi,o)(M2,®  -  M2,o) 
(1  ~  Mo)  .  0-1,0  CT2,0  O’it0O’2,0 


Minimizing  the  above  deflection  with  respect  to  (<jfi(*)i02(‘))  we  obtain  the  optimal  nonlinearities 
of  the  sensors. 

In  this  section,  the  situation  in  which  the  univariate  and  the  second-order  joint  densities 
of  the  two  observation  sequences  ( k  ~  1,2)  are  identical,  i.e.,  fi(x)  =  f2{x)  =  f(x) 

and  fi\x,y)  =  f\ 2\x,y)  —  f^\x,y),  for  j  =  1,2,  is  of  particular  interest.  For  these 

symmetric  densities,  the  optimal  nonlinearities  will  have  the  same  form,  i.e.,  (<7i(-)>M2(-))  satisfy 
gi(-)  =  g2{ •)  =  </(•)•  Under  the  these  symmetric  (identical)  conditions,  the  means  and  the 
variances  of  the  nonlinearities  satisfy 


Me  —  Mi,e  =  M2,® 


as  well  as 


and,  as  6  — ►  0,  (57)  takes  the  form 


oo  =  o-i,o  =  02,0 


rm  2(/.,-„o)2  2t>W  .... 

In  this  situation  of  symmetric  conditions,  we  can  prove  the  necessary  and  sufficient  conditions 


for  the  optimal  nonlinearity  (see  [15]).  For  the  general  situation  of  asymmetry,  however,  a 
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numerical  test  has  to  be  conducted  by  using  the  continuous-time  optimization  techniques.  The 
necessary  and  sufficient  conditions  for  the  minimization  of  (60)  with  respect  to  g(-)  has  been 
accomplished  for  m-dependent  noise,  </>-mixing  or  p-mixing  noise  in  [15]  and  will  not  be  repeated 
it  here.  The  final  form  of  linear  integral  equation,  which  determines  the  optimal  <?(•),  is 

/oo 

Kc(x,y)g(y)dy=  g(x)  (61) 

-OO 

where  Kc(x,y)  is  the  kernel  of  integration  and  has  the  form 
Kc(x,y)=  (fNW/Nw(y\x)  +  fffU)/Ni2)(v\x))/2  -  2/(2/) 


m 


with  the  notation  fN(t)  ^U)(’l-)  representing  the  conditional  density  of  xj+j  given  xj1^. 

IV. B.  Directly  Quantized  Transmission 

As  discussed  in  Subsection  III.B,  for  practical  purposes,  we  modify  the  scheme  previous 

subsection  and  adopt  the  same  structure  of  quantized  transmission  as  the  one  in  Subsection 

III.B  for  the  present  case  with  symmetric  conditions. 

By  directly  changing  the  integral  equation  (61)  to  its  discrete-form  we  obtain 

M 

~  f{xti)  ~  Kc(xt„xt})g{xt})Axt}  =  g(xti)  (63) 

i=o 

where  £<t  (i  =  0, 1,  •  •  •,  MO  are  the  breakpoints.  Define  the  vectors 

£=  [5(^o)  ■••,9{xtM)\  (64) 

L  -  l/K)/K)r-./'(%)]  (65) 

and  the  matrix 

G  =  [Gjj]  (66) 
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where 


Gij  =  f(xti)[Kc(xu,xt])Axtj  +  6(U  -  tj )]  (67) 

and  let  Xti  ( i  =  0, 1,  •  •  •,  M)  be  chosen,  as  that  the  matrix  G,  is  nonsingular.  Then  we  can  solve 
(63)  for  £  as 

gT  =  G- lfT.  (68) 

From  the  discrete-form  g(xt{ )  (t  =  0,1,*  •  •,  M)  the  quantizer  for  both  sensors  takes  the  following 
form  for  any  observation 

if  x  <  x\^ 

S(x)=i  [g(xti)  + g{xu+1)]/2  if  x  e  (*i,fc), *«*+,]  (69) 

'{x>  XtM-l 

for  i  =  1,2 ,  •  •  -  ,M  -  2.  This  will  be  transmitted  to  the  fusion  center  n  times  to  form  the  test 
statistic  Tn,k  —  g(x\^)  (&  =  1,2),  which  will  be  combined  as  in  (54)  for  computing  the 

likelihood  ratio  function. 

IV.C.  Optimally  Quantized  Transmission 

We  consider  the  scheme  of  in  Subsection  III.C  for  the  present  case.  Under  the  symmetric 
densities,  we  have  fiitg  =  fi2,o  =  fie,  for  all  6  >  0;  =  02,0  =  #o>  and  the  optimal  quantizers 

for  the  two  sensors  have  the  same  form  because  of  the  symmetricity,  i.e.,  Q\  —  Q%  =  Q .  Denote 
by  t  =  •  •  *»Uv/]  the  breakpoints  and  by  u  =  [«i,  «2,  *  •  •,  «m]  the  levels  of  Q.  From  the 

symmetric  second- joint  densities  we  have 

JP0-{-X'P)  G  (<i-i,*.*],-Xji\  €  (*i_i,*i]> 

=  P0{X[2)  e  (*,■_!, €  (<;_!,*,]} 

=  Po{Xt  6  (tj-i,t,], Xj+i  6  (tj-i,tj]} 
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Let  p(Q )  be  the  correlation  coefficient  of  fn<1  and  fn>2  and  /  the  joint  density  of  n[1]  and  Jvj2\; 

then 


ftQ)rt(Q)  =  E{Q(x[1))Q{x[2))]+2^E{Q{x[1))Q(X]%)}  -  (2m  +  l)(£[Q(xi1})])2 

MM  rt  rti  m  M  M  .  . 

=  EEw  /  ’  /  />,2/)^2/+2EEEti*Ui •  Po^i  e  .xjVi  e 

t_l  {—I  Jti— 1  Jtl—l  j= 1  1=1  i=l 

M 

-(2m  + 1)[^«/  /  /(:r)da;]2,  (70) 

/=i  Jt‘~1 

and  the  efficacy  for  the  design  criterion,  which  drives  from  the  deflection  in  the  fusion  center, 


has  the  form 


E(C)  -  2  (n^'o(^2 

MCJ)-J(l  +  p(Q))nag(Qy 


(71) 


Define  the  matrix 


Pc  =  tftj 


(72) 


with 

Pij  =  f  f  f(x,y)dxdy  +  2  £><*!"  6  6  (*,-„!,]}.  (73) 

Jti-1  Jtl _J  j  =  l 

Let  Fi  =  F2  =  F,  P\  =  P2  =  P  (where  F*  and  Pk  are  defined  in  Subsection  III.C)  and 
p!  =  R2  =  F.  Moreover,  let  us  assume  that  the  ti  (/  =  0, 1,’**,M)  are  chosen,  such  that  the 
matrix  F  +  P  +  P  —  (4m  +  2 )R  is  positive  definite.  Following  steps  similar  to  the  ones  in  [9]  we 
compute  the  optimal  quantization  levels  used  by  sensor  k  for  fixed  breakpoints 

«=  -(F  +  P  +  Pc-(4m  +  2)R)-1Al  (74) 

where  Af  =  f{ti)  -  f(U- i)J  =  1,2,---,M.  The  efficacy  employing  the  optimal  levels  at  the 
fusion  center  is  a  function  of  the  breakpoints  and  has  the  form 

E  =  AfT(F  +  P  +  Pc-  (4m  +  2)R)~1Af_  (75) 
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where  we  omit  the  arguments  i,  (/  =  0,1  Again,  the  objective  function  (75)  can  be 

optimized  with  respect  to  the  breakpoints  by  using  the  numerical  optimization  techniques. 
IV.D.  Optimally  Quantized  Transmission  with  Nonlinear  Preprocessors 

Here  we  consider  the  counterpart  of  Subsection  III.D  for  the  present  case.  Under  the  sym¬ 
metric  densities  of  Subsection  IV.A,  the  optimal  nonlinearities  ( gi(-),g2( •))  satisfy  fin(-)  = 
g2(.)  =  g(.y,  thus  =  P2,e  =  for  a11  #  >  °,  and  ^1,0  =  °2,o  =  <r0.  The  quantizers 

(9i(.),92(-))  satisfy  qi(-)  =  q2(-)  =  q(‘)\  thus  Mi =  M2,e{9,q)  =  Me(g,q),  for  6  >  0,  and 
Vio(g,q)  =  V2,o(g,q)  =  V0{g,q).  The  design  criterion  is  from  the  maximization  of  the  efficacy 

with  respect  to  ( g,q ),  i.e., 

X  2  _ (fqiVnvox  +  Wo)f'G(x)dx )2 _  (76) 

mg^a0’  ’  (1  +  />,)[JY(Vw  +  nti0)fG{x)dx  -  {J  q(y/Ea0x  +  n^)fG(x)dxf\ 

where  pq  is  the  correlation  coefficient  of  q(Tnf i)  and  q{Tn ,2). 

By  using  techniques  similar  to  the  ones  of  Subsection  III.D,  the  above  complex  optimization 
problem  can  be  separated  in  two  decoupled  optimization  problems.  The  one  with  respect  to  g 
has  the  solving  procedure  discussed  in  Subsection  IV.A,  the  one  with  respect  to  q  was  solved  in 
Subsection  IV.C  for  the  sample  size  n  -  1. 
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V.  NUMERICAL  RESULTS 


To  illustrate  the  analysis  and  the  performance  of  the  quantization-fusion  schemes  devised 
in  the  previous  sections,  we  consider  the  weak  signal  in  Cauchy  noise  which  is  symmetric  with 
respect  to  its  medium  for  the  specific  case  of  two-sensor  data  fusion.  We  consider  p- mixing  noise 
truncated  to  m  which  is  large  enough  for  the  simulation  and  the  sample  size  is  10s.  Without 
loss  of  generality,  the  the  weak  signal  8  is  assumed  to  be  0.06  for  convenience.  The  univariate 
density  of  the  Cauchy  noise  in  our  numerical  examples  is  given  by 

1 


h(x^)  = 


r,  *  =  l,2 


7r[l  +  (x(k)  -  7fc)2]’ 
where  — oo  <  7*.  <  00  is  the  median  (assumed  known). 

Although  the  second-order  joint  density  of  a  Cauchy  noise  is  difficult  to  characterize  directly, 
it  can  be  calculated  from  the  second-order  joint  density  of  a  Gaussian  process  by  a  nonlinear 
transformation  (see  [6]).  Let 

exp[-(x(fc)-7jt)2/2] 


fc(x^)  = 


x/27 


and 


/o(»(*Wfc))  =  -(ij-  exp{- — \^[{x{k)-ik)2Hy(k)-ik)2-^pG\x[k)~ik)(y(k)-ik)}} 

Ml  -Pa)  2(1  —  Pq) 

be  the  univariate  and  the  second-order  joint  densities  of  the  underlying  Gaussian  process,  where 
Pq^  is  the  correlation  coefficient.  Then  the  nonlinear  transformation  mentioned  above  has  the 
form  T(x)  =  tan[7r  •  er/(x/\/2)/2]  (see  [9]).  In  addition,  we  assume  that  the  underlining  Gaus¬ 
sian  process  for  the  corresponding  Cauchy  noise  of  each  sensor  is  characterized  the  following 
autoregressive  model  with  the  correlation  coefficients  -1  <  pk  <  1  for  k  =  1,2 


N(k)  _  VW 

Nlk)  =  PkN'k\  +  ^1  -  p\vlk)-  1  >  1 
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where  both  V^k\  for  i  —  1,2,  •  •  -  ,n;  k  =  1,2,  are  sequences  of  i.i.d.  Gaussian  random  variables 
and  have  standard  Gaussian  densities.  In  the  following  examples,  for  the  cases  of  dependence 
across  time  only,  for  i  =  1, 2,  •  •  • ,  n  are  generated  independently  of  each  other;  for  the 

cases  of  dependence  across  time  and  sensors,  they  are  generated  dependently  as  follows 

V}2)  =  pcVt{l)  +  yJl-plWi 

where  and  W{  are  two  independent  i.i.d.  Gaussian  processes  and  —  1  <  pc  <  1  is  another 
correlation  parameter. 

In  the  following  Tables  for  the  quantization  levels  and  breakpoints,  only  the  right-half  (pos¬ 
itive  if  the  medium  of  the  Cauchy  noise  is  zero)  values  of  the  breakpoints  and  the  revalent 
level  values  are  given  since  the  resulted  quantizers  (levels  and  breakpoints)  are  symmetric  with 
respect  to  the  medium  of  the  noise. 

V.A.  The  Case  of  Correlation  Across  Time  Only 

In  the  examples  for  the  case  of  dependence  across  time  only,  we  set  71  =  0.00,  72  =  0.50, 
Pi  =  0.95  and  pi  =  0.90;  the  truncated  m  is  150  for  sensor  1,  and  100  for  sensor  2. 

Example  1:  This  example  is  for  the  scheme  of  unquantized  transmission  in  Section  III.  By 
employing  the  above  values  of  parameters,  the  optimal  nonlinearities  gi,0pt(z)  and  g2,opt{x )  are 
plotted  in  Figures  1.1  and  1.2,  in  which  plots  of  the  corresponding  locally  optimal  nonlinearities, 
gk,Ud  (k  =  1,2),  obtained  by  ignoring  the  dependence  are  also  included.  In  Figure  1.3,  the 
ROCs  for  the  schemes  of  employing  these  nonlinearities  are  given,  and  the  ROC  of  single-sensor 
detection  by  using  gil0pt  is  also  included  there  (represented  by  dotted  lines). 

Example  2:  This  example  is  for  the  scheme  of  directly  quantized  transmission  in  Section 
III.  Since  the  values  of  the  Cauchy  random  variables  for  the  noise  process  are  transformed 
from  the  underlying  Gaussian  random  variables,  the  range  of  the  breakpoints  (the  interval  of 
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integration)  is  set  for  these  Gaussian  random  variables.  Thus,  we  set  (Xf*}n,  Xmlx)  =  (—5.0, 5.0) 
and  (xj£}n,Xmax)  =  (-4.5, 5.5),  which  are  enough  for  the  Gaussian  random  variables  with 
the  mean  zero  and  the  variance  one,  as  well  as  with  the  mean  0.50  and  the  same  variance 
one,  respectively.  For  the  purpose  of  good  performance  and  the  consideration  of  practice,  the 
numbers  of  quantization  levels  are  reasonably  chosen  to  be  M  =  32  for  both  sensors.  The 
resulted  quantization  levels  and  breakpoints  for  the  Cauchy  noise  are  given  in  Table  1. 
Example  3:  In  this  example,  the  values  of  the  parameters  for  each  sensor  are  the  same  as  in 
the  previous  one,  but  are  used  for  the  scheme  of  optimally  quantized  transmission  in  Section 
III.  The  optimal  quantization  levels  and  breakpoints  are  given  in  Table  2. 

Example  4:  This  example  is  for  the  scheme  of  optimally  quantized  transmission  with  pre¬ 
processors  in  Section  III.  The  optimal  nonlinearities  are  the  same  as  in  Example  1;  the  optimal 
quantizers  are  characterized  in  Table  3,  where  the  breakpoints  use  the  scale  of  standard  Gaussian 
random  variables. 

V.B.  The  Case  of  Dependence  Across  Time  and  Sensors 

In  the  examples  for  the  case  of  dependence  across  time  and  sensors,  we  set  =  72  =  7  =  0.0, 
pi  —  p2  =  p  =  0.95  and  m  =  200  for  both  sensors. 

Example  5:  This  example  is  the  for  the  same  scheme  as  in  Example  1  but  of  this  case.  The 
optimal  nonlinearity  gopt  used  by  both  sensors  are  drawn  in  Figure  5.1,  in  which  the  optimal 
nonlinearities,  gidd  and  gad,  obtained  by  ignoring  the  dependence  across  sensors  and  across  time 
and  sensors,  respectively,  are  also  included.  In  Figure  5.2,  we  draw  the  ROCs  of  employing  these 
nonlinearities;  the  ROC  of  single-sensor  detection  is  also  included  there  (represented  by  dotted 
lines). 

Example  6:  This  is  the  corresponding  example  for  the  same  scheme  as  in  Example  2  but  of 
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current  case.  The  same  parameters  for  the  underlying  Gaussian  process  (M  =  32,  7  =  0.0, 
Xmin  =  —5.0  and  Xmax  —  5.0)  as  of  sensor  1  in  Example  2  are  set.  The  optimal  quantizer  are 
given  in  Table  4. 

Example  7:  By  using  the  same  set  of  parameters  as  in  Example  6,  the  optimal  quantizer  for 
the  same  scheme  as  in  Example  3,  but  of  this  case,  are  given  in  Table  5. 

Example  8:  The  final  example  is  for  the  same  scheme  as  in  Example  4  for  this  case.  The 
optimal  nonlinearity  has  been  given  in  Example  5,  and  the  optimal  quantizer  is  given  in  Table 
6,  where  the  breakpoints  use  the  scale  of  standard  Gaussian  random  variables. 
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V.  CONCLUSIONS 


The  problem  of  memoryless  data  fusion  from  multiple  sensors  for  detection  of  a  weak  sig¬ 
nal  in  dependent  noise,  which  has  a  stationary  univariate  and  second-order  densities  for  each 
sensor,  was  formulated  and  analyzed.  The  design  criterion  for  the  structure  of  each  sensor  is 
the  deflection  which  is  from  the  Neyman  Pearson  test  of  Gaussian  data  in  the  fusion  center. 
By  considering  ideal  channel  capacities,  we  first  derived  the  optimal  memoryless  nonlinearities 
of  continuous-time  for  the  structures  of  the  sensors.  Then  to  fully  utilize  the  channel  capaci¬ 
ties  under  bandwidth  consideration  three  distinct  quantization  schemes  were  proposed  for  the 
information  transmission  between  the  sensors  and  the  fusion  center. 

In  the  first  quantization  scheme,  the  quantizers  used  in  the  sensors  were  obtained  directly 
from  the  discrete-form  of  optimal  nonlinearities,  and  were  not  optimal  quantizers.  In  the  second 
scheme,  by  replacing  the  role  of  the  nonlinearity  with  a  quantizer-type  nonlinear  function  for  each 
sensor,  we  derived  the  formula  for  the  optimal  quantizer  (optimal  breakpoints  and  quantization 
levels)  of  each  sensor.  Finally,  the  scheme  with  a  nonlinear  preprocessor  before  quantizing  for 
each  sensor  was  devised;  in  particular,  we  proved  that  the  complex  optimization  problem  with 
respect  to  the  nonlinear  preprocessor  and  the  quantizer  for  each  sensor  could  be  decoupled  in 
two  optimization  problems:  the  one  is  for  the  optimal  nonlinear  preprocessor,  and  the  other  is 
for  the  optimal  quantizer. 

In  our  work  we  assumed  that  first-  and  second-order  densities  of  the  noise  processes  are  known 
a  priori.  However,  in  most  situations  the  statistics  of  the  observations  are  not  completely  known 
beforehand  and  thus  a  robust  scheme  is  to  be  developed,  which  is  under  current  investigations. 
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sensor  1 

sensor  2 

. 

breakpoints 

levels 

breckpoints 

levels 

0.0 

0.0574 

0.5000 

0.1108 

0.4057 

0.1178 

0.9057 

0.2252 

0.9043 

0.0757 

1.4043 

0.1515 

1.6405 

-0.0051 

2.1405 

0.0052 

2.9014 

-0.0459 

3.4014 

-0.0862 

5.3253 

-0.0454 

5.8253 

-0.0989 

10.4401 

-0.0331 

10.9401 

-0.0725 

22.1622 

-0.020 

22.6621 

-0.0423 

51.2539 

-0.0102 

51.7539 

-0.0210 

129.5022 

-0.0049 

130.0022 

-0.0092 

358.0428 

-0.0020 

358.5428 

-0.0035 

1084.3265 

-0.0007 

1084.8265 

0.0019 

3600.0892 

-0.0002 

3600.5892 

-0.0003 

13112.7120 

-0.00005 

13113.2120 

-0.00008 

52426.3876 

0.000004 

52426.8876 

-0.00002 

230200.9010 

0.00003 

230201.4010 

0.000004 

Table  1:  the  quantizers  for  Example  2 


sensor  1 


sensor  2 


breakpoints 

levels 

breakpoints 

levels 

0.0 

-0.3197 

0.5000 

-0.4326 

0.0626 

-0.2297 

0.0580 

-0.3266 

0.1277 

-0.1140 

0.1352 

-0.2000 

0.2523 

-0.0170 

0.2287 

-0.0711 

0.3219 

0.0701 

0.3326 

0.4383 

0.5302 

0.1278 

0.4555 

0.1511 

0.6863 

0.1352 

0.6722 

0.2070 

0.8584 

0.1094 

0.8627 

0.1974 

1.2264 

0.0692 

1.2165 

0.1462 

1.4243 

0.0272 

1.4664 

0.0781 

2.0110 

-0.0158 

1.9771 

-0.0042 

2.6446 

-0.0406 

2.6552 

-0.0777 

3.8566 

-0.0503 

4.6643 

-0.1082 

7.0958 

-0.0433 

7.1017 

-0.0955 

13.9683 

-0.0174 

14.0342 

-0.0398 

235.7976 

0.0070 

235.3890 

1.6089 

Table  2:  the  optimal  quantizers  for  Example  3 


for  both  sensors 


breakpoints 

levels 

0.0 

0.0659 

0.1320 

0.1981 

0.2648 

0.3314 

0.3991 

0.4668 

0.5395 

0.6050 

0.6761 

0.7473 

0.8210 

0.8947 

0.9718 

1.0490 

1.1300 

1.2120 

1.2990 

1.3870 

1.4820 

1.577 

1.682 

1.788 

1.908 

2.029 

2.174 

2.319 

2.050 

2.692 

2.977 

3.263 

Table  3:  the  optimal  quantizers  for  Example  4 


for  both  sensors 


breakpoints 

levels 

0.0 

0.1152 

0.4057 

0.2735 

0.9043 

0.2720 

1.6405 

0.1759 

2.9014 

0.0937 

5.3253 

0.0446 

10.4401 

0.0188 

22.1622 

0.0079 

51.2539 

0.0034 

129.5022 

0.0020 

358.0428 

0.0017 

1084.3265 

0.0018 

3600.0892 

0.0019 

13112.7120 

0.0021 

52426.3876 

0.0023 

230220.9010 

0.0029 

Table  4:  the  quantizers  for  Example  6 


for  both  sensors 


breakpoints 

levels 

0.0 

-1.7900 

0.0672 

-1.4464 

0.1296 

-1.1168 

0.2036 

-0.7633 

0.2946 

-0.3963 

0.4075 

-0.0801 

0.5145 

0.2034 

0.6847 

0.4336 

0.8634 

0.6015 

1.4832 

0.6002 

1.8523 

0.5001 

2.7674 

0.2499 

8.5091 

0.0323 

13.9745 

-0.0378 

38,0959 

-0.0791 

1998,2669 

-0.0609 

Table  5:  the  optimal  quantizers  for  Example  7 


for  both  sensors 


breakpoints 

levels 

0.0 

-0.0686 

0.1188 

0.2874 

0.2492 

0.6231 

0.4093 

0.9254 

0.8322 

0.8917 

1.0154 

0.7492 

1.1926 

0.5852 

1.3842 

0.4230 

1.5884 

0.2862 

1.7944 

0.1850 

1.9979 

0.1153 

2.1993 

0.0695 

2.3998 

0.0404 

2.6000 

0.0227 

2.8000 

0.0123 

3.0000 

0.0042 

Table  6:  the  optimal  quantizers  for  Example  8 


Fig.  1.1  The  Nonlinearities  of  Sensor  1  in  Example 
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Fig.  1.2  The  Nonlinearities  of  Sensor  2  in  Example 
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Fig.  5.1  The  Nonlinearities  of  the  Two  Sensors  in  Example  5 
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